Scattering length for composite bosons in the BCS-BEC crossover 
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The present manuscript concerns the calculation of the boson-boson scattering length for the 
composite bosons that form as bound-fermion pairs in the strong-coupling limit of the BCS-BEC 
crossover. The material presented in this manuscript was already published as a part of a longer 
paper on the BCS-BEC crossover problem (P. Pieri and G.C. Strinati, Phys. Rev. B 61, 15370 
(2000)). Given the recent experimental advances on the formation of ultracold bosonic molecules 
from a Fermi gas of atoms by a Feshbach resonance, the calculation of the boson-boson scattering 
length has now become of particular interest. The present short version of the above paper could 
thus be helpful to the scientific community working on ultracold atomic gases. Accordingly, the 
present manuscript is intended for circulation as a preprint only. 



The recent formation of ultracold bosonic molecules 
from a Fermi gas of atoms 1 by a Feshbach resonance al- 
lows for an experimental check of theoretical calculations 
for physical quantities within the BCS-BEC crossover. In 
particular, in Ref. 2 the relation between the composite- 
boson scattering length as and the fermionic scattering 
length ap was calculated in the strong-coupling limit 
of the BCS-BEC crossover. The summation therein 
of all bosonic T-matrix diagrams has led to the result 
as = 0.75ai?- This result corrects the value ob = 2ap 
obtained within the Born approximation for the effective 
residual bosonic interaction 3-5 . The result as — 0-75ap 
could be tested experimentally in the near future, by 
measuring at the same time the molecule-molecule scat- 
tering length and the fermion-fermion scattering length 
while scanning the magnetic field through the Feshbach 
resonance. 

In this manuscript, we provide a condensed version of 
the material published in Ref. 2, focusing specifically on 
the calculation of the bosonic scattering length. We hope 
that this short summary of our previous work could be 
useful to the scientific community at the present time. 



I. BUILDING BLOCKS OF THE 
DIAGRAMMATIC STRUCTURE FOR 
COMPOSITE BOSONS 

In this section, we discuss the diagrammatic structure 
that describes generically the composite bosons in terms 
of the constituent fermions. Our theory rests on a ju- 
dicious choice of the fermionic interaction, which (with- 
out loss of generality) greatly reduces the number and 
considerably simplifies the expressions of the Feynman 
diagrams to be taken into account. 



A. Regularization of the fermionic interaction 

We begin by considering the following Hamiltonian for 
interacting fermions (we set Planck h and Boltzmann ks 
constants equal to unity throughout): 
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Here, ip a {T) is the fermionic field operator with spin pro- 
jection a = (t, |), m the fermionic mass, /i the fermionic 
chemical potential, and V c s(r — r') the effective potential 
that provides the attraction between fermions. For the 
application to atomic gases, the two spin states corre- 
spond to two different hyperfine states of the fermionic 
atoms. 

To simplify the ensuing many-body diagrammatic 
structure (yet preserving the physical effects we are af- 
ter) , we adopt for V e g the form of a "contact" potential 6 



y eff (r - r') = v S(r - r') 



(1.2) 



where vq is a negative constant. With this choice, the in- 
teraction affects only fermions with opposite spins in the 
Hamiltonian (1.1) owing to Pauli principle. A suitable 
regularization of the potential (1.2) is, however, required 
to get accurate control of the many-body diagrammatic 
structure. In particular, the equation (in the center-of- 
mass frame) 
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for the fermionic scattering length ap associated with 
the potential (1.2) is ill-defined, since the integral over 
the three-dimensional wave vector k is ultraviolet diver- 
gent. The delta-function potential (1.2) is then regular- 
ized, by introducing an ultraviolet cutoff fco in the inte- 
gral of Eq. (1.3) and letting vq — > as fco — ► oo, in order 
to keep ap fixed at a finite value. The required relation 
between vq and fco is obtained directly from Eq. (1.3). 
One finds: 
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when k \a F \ S> 1. 

With the regularization (1.4) for the potential, the clas- 
sification of the many-body diagrams gets considerably 
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simplified, since only specific sub-structures of these di- 
agrams survive when the limit fco — > oo is eventually 
taken. In particular, in order to obtain a finite result 
for a given Feynman diagram, the vanishing strength w 
of the potential should be compensated by an ultraviolet 
divergence in some internal wave-vector integration. For 
the particle-particle ladder of Fig. 1, the internal wave- 
vector integration associated with each rung diverges in 
the limit fco — ► oo and compensates the vanishing of vo, 
yielding the finite result: 
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Here Q°(k) — [iu> n — k 2 /(2m) — is a bare fermionic 
Green's function (u„ being a fermionic Matsubara fre- 
quency). Note that the interaction (1.7) depends on wave 
vectors as well as on Matsubara frequencies, revealing in 
this way the composite nature of the bosons. The factor 
of 2 in the definition (1.7) corresponds to the two differ- 
ent sequences of spin labels that can be attached to the 
four fermionic Green's functions, as shown graphically in 
Fig. 2(b). 
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Here, £(k) = k 2 /(2m) — /i and q = (q, is a four- 
momentum, with wave vector q and Matsubara fre- 
quency (v integer). In an analogous way, one can 
show that in the particle-particle channel the contribu- 
tions of the vertex corrections and of the two-particle 
effective interactions other than the rung vanish for our 
choice of the potential. 

It is thus evident from these considerations that, with 
our choice of the fermionic interaction, in the strong- 
coupling limit the skeleton structure of the diagrammatic 
theory can be constructed only with the particle-particle 
ladder (1.5) plus an infinite number of interaction ver- 
tices. A careful diagrammatic analysis considered in de- 
tail in Ref. 2 then shows that: (a) Bare composite bosons 
correspond to fermionic particle-particle ladders; (b) The 
interaction among bare composite bosons is described by 
4-point, 6-point vertices, and so on, which correspond to 
the product of 4, 6, . . ., fermionic bare Green's functions 
(with one internal four-momenta integration). The cor- 
respondence rules for the bosonic Green's function and 
the 4-point vertex are shown in Fig. 2. 

In particular, in the strong-coupling limit (whereby 
(3\fi\ 3> l), 7 the particle-particle ladder Fo(g) has the 
following polar structure: 3 
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where we have used the definition (1b = 2,u + e for 
the bosonic chemical potential (eo = l/(ma 2 F ) being the 
bound-state energy of the fermionic two-body problem) . 
Note that (apart from the residue being different from 
unity) the expression (1.6) resembles a "free" boson prop- 
agator with mass 2m. The (four-point) effective two- 
boson interaction reads instead 
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FIG. 1. Particle-particle ladder. Four-momenta are indi- 
cated and spin labels are represented by up and down arrows. 
Dotted lines represent the interaction potential and full lines 
the fermionic bare single-particle Green's functions. 
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FIG. 2. Graphical correspondence (a) between the bare 
propagator for composite bosons (represented by a thick line) 
and the fermionic particle-particle ladder, and (b) between 
the effective two-boson interaction and the four-point vertex. 

In the strong-coupling limit /3|/i| » 1, a typical value of 
the two-boson effective interaction is obtained by setting 
qi = ■ ■ ■ = qi = in Eq. (1.7). One gets: 
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The factor m 2 aF / '(8tt) in Eq. (1.8) reflects the difference 
between the true bosonic propagator and the particle- 
particle ladder in the strong-coupling limit [cf. also 
Eq. (1.6)]. Owing to this difference, u 2 (0) given by 
Eq. (1.9) (and not u 2 (0) given by Eq. (1.8)) has to be 
identified with the boson-boson interaction at zero four- 
momenta. Note that u 2 (0) is positive in the strong- 
coupling limit, thus ensuring the stability of the bosonic 
system. 

Recall further that the scattering length a^ orn within 
the Born approximation, obtained for a pair of true 
bosons (each of mass 2m) mutually interacting via a two- 
body potential with Fourier transform u 2 (0) at zero wave 
vector, is given by a®° rn = 2mu 2 (0) / (4ir) . Equation (1.9) 
then yields the following relation between the bosonic and 
fermionic scattering lengths within the Born approxima- 
tion: 
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The result (1.10) was also obtained in Ref. 3 within the 
fermionic self-consistent T-matrix approximation (which 
corresponds to the bosonic Hartree-Fock approximation 
in the strong-coupling limit). In that reference, the result 
(1.10) was erroneously regarded to be the full value of the 
scattering length a B for a "dilute" system of composite 
bosons. We will, in fact, show in the next section that the 
result (1.10) differs from a B when all bosonic T-matrix 
diagrams are properly taken into account. 

Besides the four-point vertex (1.7), the composite na- 
ture of the bosons produces (an infinite set of) additional 
vertices. One can show that all interaction vertices can 
be neglected in comparison with the four-point vertex in 
the strong-coupling limit. 2 In this limit, one can thus con- 
struct all diagrams representing the two-particle Green's 
function in the particle-particle channel in terms only 
of the bare ladder and of the four-point interaction ver- 
tex. This is precisely what one would expect on physi- 
cal grounds, since the interactions involving more than 
two bodies become progressively less effective as the com- 
posite bosons overlap less when approaching the strong- 
coupling limit. 



II. NUMERICAL RESULTS FOR THE 
SCATTERING LENGTH OF COMPOSITE 
BOSONS 

In three dimensions the scattering length a character- 
izes the low-energy collisions for the scattering from an 



ordinary potential. For the mutual scattering of two par- 
ticles (each of mass M), a can be expressed by the re- 
lation t(0) = 4ira/M in terms of the ordinary T-matrix 
t(0) in the limit of vanishing wave vector. In particular, 
within the Born approximation t(0) is replaced by the 
Fourier trasform u(0) of the interparticle potential for 
vanishing wave vector. 

For composite bosons, the T-matrix t B (qi,q 2 , 93, 94) is 
defined by the following integral equation (cf. Fig. 3): 

*b(91,92,93,94) = U2(9l,92,93,94) 

~7jy «2(9i,92,95,9i+ 92 - 95)r (95)r (9i + 92 - 95) 

*t B (qi + 92 - 95,95,93,94) (2.1) 

where u 2 is proportional to the effective two-boson inter- 
action of Eq. (1.7): 
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xg°(-k)g°(k + q 2 )g°(-k + qi - q4 )g°(k + Q4 ) . (2.2) 



FIG. 3. T-matrix diagrams for the composite bosons. 
Heavy lines and squares are denned as in Fig. 2. 

Similarly to the problem of point-like bosons, we de- 
fine the scattering length as for the composite bosons 
(each of mass 2m) in the strong-coupling limit and for 
vanishing density, by setting t B (0) = Anas/ (2m). Here 
t B (0) = (8Tr/(m 2 a F )) 2 t B (0) and ?s(0) = i B {0, 0,0,0). 
The rescaling between ts(0) and fs(0) is due to the 
difference between the true bosonic propagator and the 
particle-particle ladder in the strong-coupling limit. The 
same rescaling is consistently used when defining the 
boson-boson interaction [cf. Eqs. (1.8) and (1.7)]. 

To the lowest order in the effective interaction between 
the composite bosons, we can replace ts(0) by u 2 (0) and 
write M2(0) = 47ra^ orn /(2m), within the Born approxi- 
mation. Comparison with Eq. (1.9) yields then the value 
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2qf, as anticipated by Eq. (1.10). 



In order to obtain the exact value of is(0) (and hence 
of the scattering length as), it is convenient to determine 
first the auxiliary quantity is (9, —9, 0, 0) by solving the 
following closed-form equation 



*b(9, -9, 0, 0) = u 2 (q, -9, 0, 0) 



(2.3) 



~W E u 2 (Q,-Q,q',-Q')ro(q')To(-q')i B (q',-q',0,0) , 
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obtained from Eq. (2.1) by setting 91 = —q 2 = q and 
q 3 = <74 = 0. This integral equation can be solved by 
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standard numerical methods, e.g., by reverting it to the 
solution of a system of coupled linear equations. Note 
that, since we are here interested in the calculation of a 
two-boson quantity the zero-density limit has to be taken 
in Eq. (2.3). This implies that \ib = (or cquivalently 
2[i = e ) and T = 0. 7 The discrete bosonic frequency Q„ 
becomes a continuous variable accordingly. 

Numerical calculation of Eq. (2.3) requires us to in- 
troduce a finite-size mesh for the variables (|q|, fi) and 
(|q'|, Q'), with the angular integral over q' affecting only 
the function U2(q, —q, q', — ?')■ Equation (2.3) is thus re- 
duced to a set of coupled equations for the unknowns 
ts(|q|, O; |q|, —ft; 0; 0), which were solved by the Newton- 
Ralphson algorithm with a linear interpolation for the 
integrals over |q'| and fi'. In this way, we are led to the 
result 

a B =0.75a F (2.4) 

within an estimated 5% numerical accuracy. 

To summarize, we have shown that, in the strong- 
coupling limit of the fermionic attraction, the value = 
0.75 a_F is obtained by a correct summation of the bosonic 
T-matrix diagrams for the composite bosons which form 
as bound-fermion pairs. This theoretical prediction could 
be tested against the experimental data with ultracold 
Fermi atoms in a trap, when bosonic molecules are ob- 
tained by a Feshbach resonance. 



[Phys. Rev. Lett. 80, 3177 (1998)] where it has been shown 
that, in the strong-coupling limit, there exists a competition 
between pair and quartet condensation in a Fermi liquid 
with finite-range attraction. It is then clear that adopting a 
finite-range fermionic potential merely makes the BCS-BE 
crossover more involved, and it is thus not relevant for the 
present purposes. 

A precise definition of what is meant by strong coupling is 
as follows. In this limit, the binding energy eo is much larger 
than the density p and temperature T. Under this condition, 
the fermionic chemical potential approaches the value — eo/2 
and /3\fj,\ > 1. 
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